Dispersion relation of the nonlinear Klein-Gordon equation through a variational 

method 



We derive approximate expressions for the dispersion relation of the nonlinear Klein-Gordon 
equation in the case of strong nonlinearities using a method based on the Linear Delta Expansion. 
All the results obtained in this article are fully analytical, never involve the use of special functions, 
and can be used to obtain systematic approximations to the exact results to any desired degree of 
accuracy. We compare our findings with similar results in the literature and show that our approach 
leads to better and simpler results. 



The nonlinear Klein-Gordon equation describes a variety of physical phenomena such as dislocations, 
ferroelectric and ferromagnetic domain walls, DNA dynamics and Josephson junctions. The simple 
sinusoidal solutions to the linear wave equation, which provide a dispersion relation independent of 
the amplitude, are lost when nonlinear terms are considered. As a matter of fact the exact solutions 
cannot be expressed in a simple form in terms of their linear counterparts, although they may still 
be oscillatory. Moreover the dispersion relation obtained in this case turns out to depend upon the 
amplitude. The solution of the nonlinear wave equation poses an interesting challenge, expecially 
in the case of strong nonlinearities, where perturbation theory by itself is not applicable: indeed in 
such cases the perturbative series does not converge and no sensible information can be extracted 
directly from it. Here we present a variational method based on the linear delta expansion to find 
fully analytical approximate dispersion relations for the nonlinear Klein-Gordon and the Sine-Gordon 
equations for weak and strong nonlinearities. Our method can be easily generalized to other cases and 
provides a systematic way to achieve the desired degree of accuracy. The solutions obtained in this 
paper are fully analytical and never involve the use of special functions. 



In this article we study the problem of describing the propagation of traveling waves obeying the nonlinear Klein- 
Gordon and the Sine-Gordon equations. Under certain conditions the effect of the nonlinearity is to preserve the 
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oscillatory behavior of the solutions and, at the same time, modify the dispersion relation for the traveling waves 
which turns out to depend on the amplitude of oscillation. As a matter of fact we are considering conservative 
systems, for which the dynamics can be mapped to the nonlinear oscillation of a point mass in a one-dimensional 
potential. The main goal of this article is to explore the effects of the nonlinearity on the solutions, providing 
simple and efficient approximations. Although for weak nonlinearities, this task can be accomplished by applying 
perturbative methods (corresponding to performing an expansion in a small parameter which governs the strength 
of the nonlinearity itself), the situation is more complicated in presence of strong nonlinearities. In such a regime 
perturbation theory cannot be applied, since the perturbative series do not converge. 

Such a problem was studied in Ref. |l| , where nonperturbative formulas for the dispersion relations of the traveling 
wave in the Klein-Gordon and the Sine-Gordon equations were derived. The formulas obtained by Lim et al. provide 
an accurate approximation to the exact results even when the nonlinearity is very strong. 

In this article we consider the same problems of Ref. lj and apply to them an approach which has been developed 
recently. QQI^IQQ Our approach is fully nonperturbative in the sense that it does not correspond to a polynomial in 
the nonlinear driving parameter and, when applied to a given order, allows us to obtain analytical expressions for the 
dispersion relations, which never involve special functions, to any desired level of accuracy. It is worth mentioning that 
in the case of weak nonlinearities, an expansion of the nonperturbative results in powers of the nonlinear parameter 
is sufficient to recover the perturbative results. 

Let us briefly describe the problem that we are interested in. We consider the nonlinear Klein-Gordon equation 

u tt - u xx + V'(u) = , (1) 

where V'(u) is a function of u, which we will assume to be odd, and the prime is the derivative with respect to u. To 
determine the periodic traveling wave, we set 



u(0), 6 = kx-Lut. (2) 



After substituting into Eq.Q we find 



n 2 n + v'(u) = o , (3) 

where f2 2 = (u 2 — k 2 ) and ii = du/d9. u(9) is periodic with period 2ir and fulfills the boundary conditions 

u(0) = A , u(0) = , (4) 

with A being the amplitude of the traveling wave. The solution of Eq.J2J) with the previous boundary conditions 
oscillates between — A and A. By integrating Eq.J3|) and taking into account Eq.(@J we obtain 

^n 2 ii 2 + V(u) = V(A) . (5) 

Considering tt 2 > we observe that 

Q = 1 (6) 

V2f A [V(A)-V(u)}-V*du 

gives the exact expression for the dispersion relation of the nonlinear Klein-Gordon equation. We neglect the case 
Q, 2 < since there is no traveling wave for this configuration. 
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This article is organized as follows: in section^ we describe the variational nonperturbative approach and apply 
it to derive approximate analytical formulas for the nonlinear Klein-Gordon equation; in section IIIII we apply our 
method to two further nonlinear equations; finally in section Hvl we draw our conclusions. 



An exact solution of Eq. can be accomplished in a limited number of cases, depending on the form of the 
potential V(u). However, when the nonlinearities due to the potential V(u) are small, it is still possible to find 
useful approximations using perturbation theory. The focus of this section will be on the opposite situation, when 
the nonlinearities are not small and a perturbative expansion is not useful. In such a case one needs to resort to 
nonperturbative methods, capable of providing the solution even in the presence of strong nonlinearities. One of such 
methods, which we will use in the present article, is the linear delta expansion (LDE) 

The LDE is a powerful technique that has been applied to difficult problems arising in different branches of physics 
like field theory, classical, quantum and statistical mechanics. The idea behind the LDE is to interpolate a given 
problem V g with a solvable one V s , which depends on one or more arbitrary parameters A. In symbolic form 
V = V S (X) + S(V g — V S (X)). 5 is just a bookkeeping parameter such that for 5 = 1 we recover the original problem, 
and for 5 — ► we can perform a perturbative expansion of the solutions of V in 8. The perturbative solution obtained 
in this way to a finite order shows an artificial dependence upon the arbitrary parameter, A, and would cancel if the 
calculation were carried out to all orders. As such we must regard such dependence as unnatural; in order to minimize 
the spurious effects of A we then require that any observable 0, calculated to a finite order, be locally independent 
on A, i.e. that 



This condition is known as the "Principle of Minimal Sensitivity" (PMS) [llj. We call Xpms the solution to this 
equation. (In the case where the PMS equation has multiple solutions, the solution with smallest second derivative is 
chosen.) We emphasize that the results that we obtain by applying this method do not correspond to a polynomial 
in the parameters of the model as in the case of perturbative methods. 

The procedure that we have illustrated is quite general and it will be possible to implement it in different ways 
depending on the problem that is being considered. In Refs. the LDE was used in conjunction with the 

Lindstedt-Poincare technique to solve the corresponding equations of motion. Our approach here is to apply the LDE 
directly to the integral of eq. JBJ) as in Refs. Q and [(J. We will consider the potential 
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V(u) = Y + — 



(8) 



The dispersion relation in this case can be obtained using Eq. © as 



n = 



ix\J 1 — [lA 1 



(9) 



2/ ( j r (l -to sin 2 0)-V2d0 ' 



with to = 



2(i+M a r 
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We consider the following approach to obtain the dispersion relation of a periodic traveling wave. This comes from 
the equation for the period of oscillations, 



T = 



zdu 



(10) 



-A yjE - V(U) 

where the total energy E is conserved and ±A are the classical turning points. 

In the spirit of the LDE we interpolate the nonlinear potential V(u) with a solvable potential Vq(u) and define 
the interpolated potential Vg(u) — Vo(u) + 8(V(u) — Vo(u)). Notice that for 8 = 1, V§{u) = V(u) is just the original 
potential, whereas for 5 = it reduces to Vo (u) . Hence we can write Eq. ifTH]) as 

" +A V2 du 



T s = 



where 



A(«) = 



a ^E -V (u) ^l + 5A(u) 
E- E -V(u) + V (u) 



E -V (u) 

Obviously E = V(A) and E a = V (A). 

We treat the term proportional to 5 as a perturbation and expand in powers of 6. This allows us to write 



(11) 



(12) 



(2n- 1)!! 
n!2^ 



(-l) n ^ 



a y/Eo-V (u)' 



zdu 



(13) 



Observe that the integrals in each order of Eq. i|13|) have integrable singularities at the turning points because 
A(±A) is finite. Assume that |A(u)| < Ao < 1 for u G [—A, A], which happens if A, the arbitrary variational 
parameter, is chosen appropiately. Then, the series l|13|) converges uniformly for \6\ < 1/Ao, which includes the case 
6 = 1. 

For the potential given in Eq. JSJ we can choose Vq(u) = 1 ~ t 2 A u 2 as the interpolating potential and hence we have 



A(u) = 



(14) 



The parameter A should be chosen to be A > y ! ^ A —^Jl + which guarantees the uniform convergence of Eci. (|13|l . 



It is straightforward to check that at first order, 
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The PMS © with O = T yields 



The period is found to be 



^PMS 



1 + A 2 



-LiA 2 



i 21 
2 



3fiA 



TpMS 
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Correspondingly, 



n 



LDE(l) 



1 + tM 2 
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(15) 



(16) 
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In Ref. 0] it was found that with this value of Xpms all the remaining terms of odd order in Eq. I|13|) vanish. Hence, 
retaining only nonvanishing contributions, the expression for the period at order N is 



yfl + 3[lA 2 

and, correspondingly, 




2vr 
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At second order we have 



^4 + 3A 2 fi 
3A 4 /i 2 (1024 + A 2 /i(1536 + 61LA 2 /i)) 
1024(4 + 3A 2 n) 4 

At third order, the dispersion relation is given by 



_ V •* Z _ ^ 



+ 

* l LDE(i) — — ; 7 5 prr - • l^J 

3AV I 385 ,4V + 560.4V (4 + 3 A 2 /*) + 1024 (4 + 3 AV 

1 + " « 

16384 (4 + 3A 2 /i) 6 

We will compare the results obtained using our method, Eqs. (|18ll . 1|21|) . and lll'L'1) with the results obtained in 
Ref. 1], where the same problem has been solved using the harmonic balance technique in combination with the 
linearization of the nonlinear Klein-Gordon equation. The findings of Ref. [l| at first order, their expression for the 
dispersion relation coincides with our Ea. H18|) . whereas at the second order they find 



/40 + 31M 2 + ^1024 + 1472^ 2 + 42 VA 4 

Si£im(2) = I/ ^ ■ (23) 

In the left-hand panel of Fig. we make a comparison of the ratios of the dispersion relations obtained from 
Eqs. l|18f) . and H21fl - I|23|) to the exact dispersion relations for [iA 2 < 0, and in the right-hand panel of Fig.^we display 
the relative error 



A = log 10 



(24) 



^exact 

for ii A 2 > 0. We can appreciate that our variational method at second order provides a smaller error than the method 
of Ref. [i| applied to the same order. The error is further reduced by using the LDE to the third order and can be 
then systematically reduced using the general formula (|19fl . 

III. FURTHER EXAMPLES 
A. Sine-Gordon model 

We now consider the Sine-Gordon model, which is governed by the potential 



V(u) — — COSM 



(25) 



6 





|xA 



FIG. 1: (Left) Ratio of the dispersion relation from Eqs. 11811 and 12H - 12lj|l exactly for /j,A 2 < and (right) relative error A 
[see eq. H24H 1 of the dispersion relations for fiA 2 S> 0. 



and which allows us to write the nonlinear Klein-Gordon equation, also known as the Sine-Gordon equation as 

n 2 il + smu = Q. (26) 
The exact dispersion relation in this case can be obtained from 

n = — (27) 

2/; /2 (l-m 2 sin 2 t)-V 2 dt 

with m = sin(^4/2). Observe that in this case 

,7r/2 

T = 4 / {1 -m 2 sin 2 ty 1/2 dt = 4K(m 2 ) , (28) 
J o 

with K(m) being the elliptic integral of the first kind. We take advantage of this fact and make use of the nonpertur- 
bative series for the elliptic integral which was derived using the LDE technique 12J. At order N, setting A = — m/2 
and 5 = 1, it is given by the expression: 

Kl ^-^^ r(j- + l/2) (-m) k 

K N[ m,A)- 2 j\?{k-j)\T{l/2-k) 2*-i (i_™)*+i/a • [Zy) 

This expression provides a nonperturbative series for the elliptic integral of the first kind since it does not correspond 
to a simple polynomial in m. 

To further improve this series we can use the Landen transformation 



1 + \ (1 + y/m) 

and the inverse relation 



KM = i ~~~i = K { ) (30) 



m 



Notice that f(m) — ^^JLy maps a value < m < 1 into a new value m! = /(to) > m. The inverse transformation 
/ _1 (m) = - 2+2 ^ n -/ n+m ' > maps a value m into a smaller one. Using this transformation we obtain more accurate 
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approximations for the elliptic integrals. For example, at order 1 we find 
and, correspondingly 



K LDE (i) (m) = = — (32) 



^LDE(l) = ^JcOs(A)+12 



<os,_ 



3 . (33) 



At second order we find 



16 cos 2 (4) (3 + 12 cos (4) + cos(^l)) 2 ^2 + 2 cos (f) sec 4 (Jj 
" Li,£ ' lJ! ~~ 2713 + 2520 cos (4) + 2580 cos (I) + 360 cos (^) + 19 cos(2^) '' >4 ' 

It is noticeable that £Ilde(3) — ^lde(2)- In fact, for the following consecutive orders, the same statement holds, 
i. e., £I>lde(S) = &lde(4), ^lde (7) = ^lde (6) an d so on. The same pattern of equal value of the observables for 
consecutive orders of approximation was found in Ref. 3] for the Duffing potential at large n. 
For comparison, Lim et al. [l| have found the dispersion relation to be given at first order as 



V-l M>=\F^, (35) 



and at second order as 



where 



and 



«£<m(2) = \/g(A) + y/g^A)-h{A) , (36) 



(6 - b 2 - bi + b 6 )A + 18oi + 2a 3 
9{A) = 38A 

h(A) = "l^-^-^ + fce) (3?) 



i 1 = 2J 1 (A), a 3 = -2J 3 (A), 6 24 = 2(-l) i J 2i (A), z = 0,l,2,3 (38) 



J n (A) being the Bessel function of the first kind. 

In the left-hand panel of Fig. [21 we display the ratio of the dispersion relations from Eqs. (|33|) - (|36|l to the exact and 
on the right-hand panel the corresponding relative errors. From the graphs we see that the LDE curves calculated to 
second order display much smaller errors than the curves obtained with the method of Lim et al. even close to A = it. 
A second observation is that our formulas can be systematically improved simply by going to a higher order and that 
they do not involve any special function, as in the case of Eq. H35JI . 



B. Pure quartic potential 



Our final example is the Klein-Gordon equation in a pure quartic potential 

V(u) = £ , 



(39) 
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A (degrees) A(degrees) 



FIG. 2: (Left) Ratio of the dispersion relation from eqs. ffl33P - 1361 1 exactly and (right) relative error A [see eq. 1241 1 of the 
dispersion relations. 



which leads to the equation of motion 

u + u 3 = . (40) 

This is a particular case of the first example where the contribution of the quadratic term in the potential JHJ) 
is neglected. As such, the corresponding dispersion relation can be derived from the expression of the period of 
oscillations, Eq. I|19l) . since the quadratic term contributes with the 4 in the square root in the front of the double 
sum and in the argument in the sum, and is simply given by 

and correspondingly Hlde(N) can be obtained as in Eq. (|20[1 . 
Results for the first three orders are the following : 

13824V3A 

^lde (2) - 28259 - (43) 
_ 1990656^3^4 

n ^(3) = 4069681 ' (44) 
Lim et al. found, at first and second order of approximation, respectively, 

x/3 



1 



n Li m{2) = —\/62 + 2V42lA. (45) 

In the left-hand panel of Fig. |21 we display the ratio of the approximate to the exact dispersion relation and in 
the right-hand panel the relative error from our findings at first, second and third order and those of Ref. |l| given 
previously. At first order, our findings perform just as the second order of Lim et al. [J, and at second and third 
orders, the performance of the variational results is excellent. 
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FIG. 3: (Left) Ratio of the dispersion relation from Eqs. 142H - 144|I . and both expressions in Eq. (I45H exactly and (right) 
relative error A (see eq. 1241 ') of the dispersion relations. 

IV. CONCLUSIONS 

We have derived analytical expressions for the dispersion relations of the nonlinear Klein-Gordon equation for 
different potentials by means of the Linear Delta Expansion. This technique is implemented by computing the period 
of oscillations in the given potential. In the particular example of the Sine-Gordon potential, where the dispersion 
relation is given in terms of elliptic integrals, we have implemented the LDE to compute such integral and, by means 
of the Landen transformation, we have obtained an improved series for the elliptic integral. We have observed that 
the expression obtained by using the first few terms in this series performs remarkably well even close to the A = n. 
We believe that our results are appealing in two respects: first in that they provide a systematic way to approximate 
the exact result with the desired accuracy, and second in that the expressions that we obtain never involve special 
functions, as in the case of Ref. l]. An aspect that needs to be underlined is that the method described in subsection 
[njprovides a convergent series representation for the dispersion relation, provided that the arbitrary parameter fulfills 
a simple condition. 
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